Growing multiplex networks 
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We propose a modelling framework for growing multiplexes where a node can belong to different 
networks. We define new measures for multiplexes and we identify a number of relevant ingredients 
for modeling their evolution such as the coupling between the different layers and the arrival time 
distribution of nodes. The topology of the multiplex changes significantly in the different cases 
under consideration, with effects of the arrival time of nodes on the degree distribution, average 
shortest paths and interdependence. 

PACS numbers: 89.75.Fb, 89.75.Hc and 89.75,-k 



Many systems are structured as networks, and their 
properties are now, after a decade of efforts, well un- 
derstood [Il-Q. However, a complex network is rarely 
isolated, and some of its nodes could be part of many 
graphs, at the same time. Examples include multimodal 
transportation networks 0, 01 , climatic systems [1] , eco- 
nomic markets energy-supply networks [loj and the 
human brain In these cases, each network has to 

be considered as being a part of a larger system in which 
a set of interdependent networks with different structure 
and function coexist, interact and coevolvc. So far net- 
work scientists have investigated these systems by look- 
ing at one type of relationship at a time, e.g. by analyz- 
ing collaboration networks and email communications as 
separate graphs. However, the structural properties of 
each of these networks and their evolution can depend in 
a non-trivial way on that of other graphs to which they 
are interconnected. For instance, two people can become 
friends after having been in touch by email for a while, 
or two researchers start a collaboration after having at- 
tended the same conference (i.e. being part of the same 
co-location network for a certain amount of time). Con- 
sequently, these systems are better represented as mul- 
tiplexes, i.e. graphs composed by M different layers in 
which the same set of N nodes can be connected to each 
other by means of links belonging to Af different classes 
or types. 

The characterization of multiplexes is still in its in- 
fancy. Some recent works have proposed suitable exten- 
sions to multi-layer graphs of classic network metrics and 
models 13,[l3l, and preliminary results show that mul- 



tiplexicity has important consequences for the dynamics 
of processes occurrin g in real systems, including rout- 

and epidemic 
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In this Letter wc propose and study a generic model 
of multiplex growth, inspired by classical preferential- 
attachment based models, in which the probability for 
a newly arrived node to establish connections to existing 
nodes in each of the layers of a multiplex is a function 
of the degree of other nodes at all layers. We define two 
new metrics to characterize the structure of multiplexes. 



namely the inter-layer joint-degree distribution and the 
node interdependence, and we study the effect of differ- 
ent attachment rules on the structure of the resulting 
network. We then focus on 2-laycr multiplexes, and we 
provide closed forms for both the degree distributions at 
each layer and the inter-layer degree-degree correlations 
corresponding to different attachment rules. We also in- 
vestigate the interplay between the attachment rule and 
the relative time-scales of edge creation at each layer, 
and we show that the presence of delays in the arrival of 
nodes at different layers can have a significant impact on 
the distribution of shortest paths and interdependence. 

More precisely, a multiplex is a set of N nodes which 
are connected to each other by means of edges belonging 
to M different classes or types. We represent each class 
of edges as a separate layer, and we assume that a node 
i of the multiplex consists of M replicas, one for each 
layer. We denote by V^"^ the set of nodes in layer a 
and by i?'"! the set of all the edges of a given type a. 
Consequently, the a*'' layer of a multiplex is represented 
by a graph G^°'\V^°'\ E^"'^), associated to the adjacency 
matrix A^"^ ~ {a!?' I- where al"' = 1 if node i and node 

j are connected by an edge of type a, whereas a["' = 
otherwise. An M-layer multiplex is fully specified by 
the vector A = [A^^^ , A^^I , • . • , A^^-'^l] , whose elements are 
the adjacency matrices of the M layers. We denote by 
~ J2j '^li' ttie degree of node i at layer a, i.e. the 
number of edges of type a of which i is an endpoint, 
and by fc^ the Af-dimensional vector of the degrees of 
the replicas of i. In general, the degrees of the replicas 
of i are distinct, and some replicas can also be isolated 
(i.e. kl"^ = for some value of a). In the following we 
consider all the edges at all layers to be undirected and 
unweighted. 

As in the case of classical 'singlex' graphs, we can char- 
acterize each layer a of a multiplex by studying the de- 
gree distribution P(fc["l), and and the joint-degree distri- 
bution P(A:1"1, 1). However, we are interested here in 
the structural properties of the multiplex as a whole, so 
we propose to quantify the correlations between the de- 
grees of replicas of the same node at two different layers a 
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and a', by constructing the inter-layer joint-degree dis- 
tributions P(fc["l, 1), or equivalently the conditional 
degree distributions P(fc[" 1 j/cl"!). In particular, we can 
look at the projection of the conditional distribution ob- 
tained by considering the average degree fc[" 1 at layer a' 
of nodes having degree fc^"' at layer a: 

fc["'l(/cH) ^ J2 /c["'lp(/c["'l|fc["l) (1) 

By plotting this quantity as a function of fc["l we can 
detect the presence and the sign of degree correlations 
between the two layers. For a multiplex with no inter- 
layer correlations between layers a and a' we expect 
A:["'l(fcM) = (fc["']) and fcM(A;M) = (/jH). If fcM(fcH) 

increases with k'^"^ we say that the degrees of the two 
layers have positive (assortative) correlations, while if 
A:["l(fc["l) is a decreasing function of fcl"! we say that 
the degrees on layer a and a' are anticorrelated (or dis- 
assortatively correlated) . 




FIG. 1: (color online) Linear attachment with syn- 
chronous arrival Panel a-b: degree distribution P{k) (left) 
and the projection k{k) of the inter-layer degree-degree cor- 
relations (right) for the {1,0,0,1} (top) and the {0,1,1,0} 
(bottom) attachment kernels (red circles are for the first layer, 
yellow squares for the second layer) . These results correspond 
to 2-layer multiplexes with A'^ — 10000 nodes, m — 3, mo — 3, 
and are averaged over 50 realizations. The observed degree 
distributions and the inter-layer degree correlations closely 
follow the theoretical curves (reported as solid black lines) 
and are relatively insensitive to the coupling matrix. 

In addition to the assortativity, we can also character- 
ize for each node i of a multiplex, its 'multiplex reach- 
ability' by computing the average distance Li from i to 
any other node of the the multiplex and comparing this 
average distance with that measured on each layer sepa- 
rately. However, the presence of more than one layer in a 
multiplex produces an increase in the number of available 
paths, so that the distance between two nodes of a mul- 
tiplex will be, in general, smaller than or at most equal 
to that measured on each layer separately. A measure 



to quantify the value added by the muitipiexicity to the 
reachability of nodes is the interdependence [l3| which 
for a node i is defined by 



where ■0^^ is the number of shortest paths between node i 
and node j which use edges lying on more than one layer, 
while <7ij is the total number of shortest paths between 
i and j in the multiplex. The interdependence of a mul- 
tiplex is computed as the average node interdependence 
A = ^/NJ2i 9,nd is normalized in [0, 1]. If A is close to 
zero, then most of the shortest paths among nodes lie on 
just one layer, while if A is close to 1 then the majority 
of the shortest paths exploit more than one layer. 

We now introduce the first model of growing multi- 
plex networks. Most of the classical growing models for 
single-layer networks assume that new nodes arrive in 
the graph one by one and that each node carries m edge 
stubs; upon arrival, the new node creates m new connec- 
tions with other existing nodes, according to a prescribed 
attachment rule. In that case, each node i has a unique 
arrival time ti. In the case of multiplexes, instead, each 
layer can exhibit a different edge-formation dynamics, 
and in general the edges of the M replicas of a new node 
are not created at the same time. For instance, a face-to- 
face interaction relationship is usually established before 
two individuals become friend, while two locations are 
usually connected by a road before a direct railway line 
between them is constructed. Consequently, we assume 
that a newly arrived node has m stubs on each layer of 
the multiplex, but the replica of a node i on layer a con- 
nects its m stubs at a different time We denote by 
ti the vector of arrival times of the replicas of node i. In 
order to make the model analytically tractable, we make 
two simplifying assumptions. The first is that there exists 
a layer a so that < i'"' \/i,\fa =/= a. This is equiva- 
lent to saying that a newly arrived node must first create 
its connections on layer a before any of its replicas can 
create connections on any other layer a ^ a. We call a 
the master layer. The second assumption is that nodes 
arrive one by one on the master layer, at equal discrete 
time intervals t = (0, 1, . . . , }. We label the nodes of a 
growing multiplex according to the ordering induced by 
their arrival on the master layer. Without loss of general- 
ity, in the following we assume that the the master layer 
is the first one, i.e. a = 1, and that the arrival times of 
the replicas of node i have the form 

iH=r(tW,eH(r)) (3) 

where T is a certain function of and of the random 
variable C'"'(t). By appropriately choosing T and £,^"^{t) 
we can model different arrival behaviours, including a) 
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FIG. 2: (color online) Linear attachment with delayed 

arrival. Blue lines correspond to /? = 1.1, red lines to 
/? = 2.0. The top sub-panels of each panel correspond to 
the {0, 1,1,0} coupling pattern, while the bottom sub-panels 
are for {1,0,0,1}. Panel a): The degree k{t) of the largest 
hub of the second layer as a function of time. k{t) scales as 
(t/toY, where 5 depends on the exponent /3 and approaches 
5 = 0.5 when /3 increases (panel b). The value of /3 tunes 
the shape of the distribution of average shortest path length 
from each node to all the others (panel c) and the distribution 
of node interdependence (panel d). When /3 increases these 
distributions become more similar to the corresponding ones 
observed for synchronous arrival (shaded grey). All plots cor- 
respond to 2-layer multiplexes with A'^ = 10000 nodes, m = 3, 
mo = 3, and the results are averaged over 50 realizations. 



simultaneous arrival (T = t'"'^'); b) power-law delayed ar- 
rival (T = + ^(t) and = t) = (^ - 1)t-'^ for 
T > 1 and l3 > 1). Upon arrival, the newborn node i 
connects to m existing nodes in the master layer, accord- 
ing to a certain attachment rule. As in the preferential 
attachment models we assume that the attachment 
probability depends on the degree of a node. However, 
in a multiplex the probability for node i to connect to 
node j on each layer a can depend not only on /sj"' but 
also on the degrees of j's replicas on the other layers 



n 



a] 



(4) 



For the sake of clarity and without loss of generality, 
we focus in the following on 2-layers multiplexes with 
a = 1,2.. 

We begin with the simplest case of linear attachment 
which is the natural extension of the Barabasi-Albert 
model 0. In this case, we consider that the probabil- 
ity for a newborn node i to connect to an existing node 
j on layer a is proportional to a linear combination of 
the degrees of j at all layers. Therefore, the attachment 
kernels can be expressed as 



FW[k,qi 



c 



[Id] c[i'2r 

[2,1] c[2'2l 



(5) 



where from now on we denote the degree fcl^l on layer 1 
by k and the degree fc'^l on layer 2 by q, i.e. A;^^! = k 
and fc[2l = q. The coefficients c^'"'*! tune the dependence 
of the attachment probability at layer r on the degrees 
of nodes at layer s. In the case of 2-layer multiplexes we 
can represent the set of coefficients C = {el'"'"!} using the 
compact notation {c^^'-^l, c^^'^l , c'^'^l, c'^'^l}. The dynam- 
ics can be easily solved in mean-field and in some specific 
cases we can also solve analytically the master-equation 
and fully characterize the degree correlations within the 
two different layers. In the following we summarize the 
master-equation solution in some particularly interesting 
cases. First of all let us consider simultaneous arrival of 
the nodes in the two layers. If we set C = {1, 0, 0, 1} then 
the attachment probability at each layer will depend only 
on the degree of the nodes in the same layer. In this case 
the degree distribution in the first layer reads 



Pik) = 



2m{m + 1) 
k{k + l){k-\-2)' 



k > 



(6) 



and the degree distribution in the second layer is identi- 
cally equal. This distribution goes as P{k) ^ k~^ with 
7 ~ 3. If we solve the master-equation for the multi- 
plex evolution we obtain the analytical expression for the 
inter-layer joint degree probability P{k,q) 



P{k,q) 



2r(2 + 2m)r(l)r(q)r(fc + q- 2m + 1) 
r(m)r(m)r(fc + q + i)T{k - m + l)T{q - m + 1) 

(7) 

The average degree k{q) at layer 1 of nodes having degree 
q at layer 2 reads: 



k{q) 



n{q + 2) 
I -\- m 



(8) 



Notice that even if the two layers grow independently, 
the simultaneous arrival introduces non-trivial inter-layer 
degree correlations. In fact, in the mean-field approach, 
the degree of a node on each layer increases over time 

as fcj"'(i) = m ' ^^^^ degrees of the two 

replicas of a node i depend, for large t, only on their ar- 
rival time. If both replicas have the same arrival time, 
i.e. tf^ = if^ then the degree of the two rephcas will 
be positively correlated. In Fig. [T] we report the degree 
distribution and the values of fc(g) for two coupling pat- 
terns, which are in good agreement with the theoretical 
curves. It is clear from the figure that in the synchronous 
arrival case the shape of the coupling matrix is actually 
not very relevant and that the value of the degree dis- 
tribution exponent and strong assortativity are robust 
features of these multiplexes. 

If we consider a power-law delayed arrival time on the 
second layer, the results are significantly different. In 
Fig. [2] we illustrate how the exponent of the delay distri- 
bution /3 affects the structure of the obtained multiplex. 
In particular, the coupling matrix has now an effect on 
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FIG. 3: (color online) Semi-linear attachment. Degree 
distributions (panel a) and inter-layer degree-degree correla- 
tions (panel b) for a multiplex with semi-linear attachment 
and different coupling patterns, respectively {1,0,0, 1} (left) 
,{0, 1, 1, 0} (middle) and {1, 0, 1, 0} (right). Red circles are for 
the first layer and yellow squares for the second layer. Panel 
c): the distribution of the average shortest path length from 
one node to all the other nodes heavily depends on the cou- 
pling pattern. Similarly, the interdependence of a node X{t) is 
always a sublinear function of the arrival time t but its shape 
depends on the coupling pattern at work (panel d). In gen- 
eral, older nodes have smaller interdependence. Panel e): The 
coupling pattern also affects the distribution of node interde- 
pendence. The smallest average interdependence is observed 
when the two layers are independent, i.e. case C = {1, 0, 0, 1} 
(yellow curve). 



various quantities. The degree distribution are still power 
laws but the exponent depends now on /3: for /3 small, 
some nodes become predominant thanks to the absence of 
other nodes and the degree distribution becomes broader. 
This effect being stronger for an off-diagonal matrix for 
which the interlayer coupling is strong. The average 
shortest path and the interdependence are also strongly 
aflFtected: the presence of more predominant 'old' hubs 
lower the average shortest path and the interlayer assor- 
tativity. 

So far, we have considered the case of two scale-free 
growing networks and it would be interesting to capture 
the possibility of having a scale- free network coupled to a 
network with peaked degree distribution. In this respect, 
we introduce the semi-linear attachment which allows to 
grow multiplexes in which the two layers have different 
topological structures. The model is defined as follows 



FW[k,qi 
F[2l[fc,g] 



C 



(9) 



where C is still a 2 x 2 matrix of coefhcients, as in the 
linear model. In this case, the degree of a node on any of 
the two layers could depend only on its degree on layer 
1 and does not ever depend on its degree on layer 2. 



If we set c^^'-'^l 



:[2.21 = 1 and c[i'2] 



.[24] 



we 



distributions of the two layers read: 
2m (m + 1) 



pW(fc) 



k{k + l){k + 2) 
1 / m 



(10) 
(11) 



1 + m \l + m ^ 
while the inter-layer joint degree distribution is equal to 



/k — m,\ 



2m 



n ) V2 + 2r?i4-fc — « 



{-If 



k — Tn-\-n 



P{k,q) = a{k) ^ 

n = ^" ' ' 

(12) 

where a{k) = r(m+i)r(k~m+i) ■ The function k{q) is 
given by 



k{q) 



2(m+ 1) 
1 + 2m 



q — m-\-l 



(13) 



Similar relations can be derived for the other coupling 
patterns. In Panel a) and b) of Fig. [3] we report the de- 
gree distribution and the value of k{q) for three different 
coupling patterns, which are in good agreeement with the 
theoretical curves. It is worth noticing that in the semi- 
linear model the coupling pattern has a dramatic impact 
on other structural properties of the multiplex such as the 
distribution of the average shortest path length from each 
node and the distribution of node interdependence. In 
particular, the interdependence is larger for older nodes, 
and grows sublinearly with time. This implies in par- 
ticular, that navigation for old nodes is easier within a 
single layer while younger nodes will have to resort to the 
different layers to reach a target. In addition, the growth 
being sublinear the performance of the system increases 
very slowly. 

In summary, we have presented a framework for mod- 
elling the growth of multiplexes which enabled us to high- 
light the relevant ingredients for these systems and which, 
we believe, will provide a solid basis for more elaborated 
models. 
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can analytically solve the master equation and the degree 
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